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Abstract 

We systematically calculate the relativistic corrections to prompt J/tf photoproduction and 
hadroproduction using the factorization formalism of nonrelativistic QCD. Specifically, we include 
the 3 S'!^ and 3 P^ color-singlet and the 3 sf^, ^q 8 ^, and 3 pf^ color-octet channels as well as the 
effects due to the mixing between the 3 s| 8 ^ and 3 d[ 8 ^ channels. We provide all the squared hard- 
scattering amplitudes in analytic form. Assuming the nonrelativistic-QCD long-distance matrix 
elements to satisfy the velocity scaling rules, we find the relativistic corrections to be appreciable, 
except in the color-singlet channel of hadroproduction. We also observe significant differences 
in the line shapes of the relativistic corrections between photoproduction and hadroproduction. 

PACS numbers: 12.38.Bx, 12.39.St, 13.85.Ni, 14.40.Pq 
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I. INTRODUCTION 


The production of heavy quarkonia, the QCD bound states of heavy-quark pairs ( QQ ), 
serves a an ideal laboratory to probe both the perturbative and nonperturbative aspects 
of QCD. The effective quantum field theory of nonrelativistic QCD (NRQCD) 1 endowed 
with the factorization formalism introduced by Bodwin, Braaten, and Lepage 2] is nowa¬ 
days the most favorable theoretical approach to study heavy-quarkonium production and 
decay. In this framework, the theoretical predictions are separated into pro cess-dependent 
short-distance coefficients (SDCs) and supposedly universal long-distance matrix elements 
(LDMEs). The SDCs may be calculated perturbatively as expansions in the strong-coupling 
constant a s , while the LDMEs are predicted to scale with definite powers of the relative 
velocity v of the heavy quarks in the quarkonium rest frame 0 . In this way, the theoretical 
calculations are organized as double expansions in a s and v. In contrast to the color-singlet 
(CS) model, in which the QQ pair must be in the CS state that shares the spin S, orbital 
angular momentum L, and total angular momentum J with the considered quarkonium, 
NRQCD accommodates all possible Fock states n = 2S+1 L^, where a = 1,8 stands for CS 
and color octet (CO), respectively. 

During the past two decades, tests of NRQCD factorization and the universality of the 
LDMEs were performed in a vast number of experimental and theoretical works. Despite 
numerous great successes, there are still some challenges in J/pj hadroproduction. As for the 
prompt yield of unpolarized J/'tp mesons, the next-to-leading-order (NLO) QCD corrections 
were calculated for all the relevant channels, including 3 5'['' 4], 3 Sf‘\ 0 , and 3 Pj 8] Q, Q 

for direct production and the feed down from ip' mesons, and 3 Pj^ 8, 9] for the feed down from 
XcJ mesons. As for the LDMEs, different ways of fitting lead to different results. By a global 
fit to the world’s J /ip data from hadroproduction, photoproduction, two-photon scattering, 


and e + e anni 
pinned down 


lilation, the three CO LDMEs relevant for direct production were successfully 


which greatly 


lo| in a way compatible with the velocity scaling rules 0], 
supported their universality. Using these LDMEs, the J/ip polarization in hadroproduction 
was predicted to be largely transversal [ill ] . By fitting to the CDF Run II measurements 
of J/ip yield and polarization for transverse momenta px > 7 GeV, the authors of Ref. [i0] 
obtained two linear combinations of the three LDMEs, which led to an almost unpolarized 
prediction for the LHC. Fitting to the prompt J/'ip yield measured for pr > 7 GeV by CDF II 
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and LHCb, including also the feed-down contributions from the XcJ and ip' mesons, a third 
set of LDMEs was obtained, which resulted in a moderately transverse Jj%b polarization at 


the LHC 


13]. All of these three LDME sets can describe well the J/ip yield at the LHC. 


Unfortunately, none of the resulting predictions for J/ip polarization is consistent with the 
latest measurements by CMS [l4] and LHCb 15]. The above results reflect the fact that 
the NRQCD prediction of prompt J/ip polarization strongly depends on the actual values of 
the LDMEs. To test the universality of the LDMEs in a meaningful way and to clarify the 
J/ijj polarization puzzle, it is useful to investigate some other effects in the determination 
of the LDMEs, such as the higher-order relativistic corrections. 

In the heavy-quarkonium system, we actually have v 2 ~ « s (2mg), which is not very small. 
In some cases, it was found that the higher-order v 2 corrections are even as important as the 
higher-order a s corrections. For example, the relativistic corrections played an important 


role in resolving both the double-charmonium m and J/V’TAnon—cc production problems at 


the B factories 


13 , 


18|]. In J /ip hadroproduction, the v 2 corrections in the CO channels 


CS 3 4 1] channel [20]. 


and :i 5'| 81 were found to be significant in the large-p^ region [1$], although they are tiny in the 

In double-quarkonium hadroproduction, the v 2 corrections also turned 




out to be significant [21], especially in the CO channels [22], In the test of the hypothesis 
A"(3872) = x'd i 11 hadroproduction, appreciable v 2 corrections were encountered in the 
channel 9]. All these observations provide a strong motivation for us to systematically study 
the v 2 corrections to the cross sections of prompt J/ip photoproduction and hadroproduction. 
This will allow us to render global fits of the contributing LDMEs more reliable and to 
deepen our understanding of their universality. While the SDCs of direct J/ip production 
immediately carry over to the feed down from the ip' mesons, the feed down from the Xcj 
mesons requires a separate calculation. 

The remainder of this paper is organized as follows. In Sec. El we explain how we calculate 
all the relevant SDCs. Our numerical results are presented in Sec. 1III1 Our conclusions are 
contained in Sec. IIV1 Our analytic results are listed in the Appendix. 


II. NRQCD FACTORIZATION FORMULA 

Invoking the Weizsacker-Williams approximation and the factorization theorem of the 
QCD parton model, the cross sections for the photoproduction or hadroproduction of the 
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hadron H = Jjif }, Xcj , ^ may be written as 10] 


<r(AB H + X) = / dx 1 dy 1 dx 2 fi/A(xi)f j / i (yi)fk/B(x 2 )a(jk ->• H + X), 

i,j,k 


( 1 ) 


where is the parton distribution function (PDF) of the parton i in the hadron A = p,p 

or the flux function of the photon i = 7 in the charged lepton A = e _ ,e + , fj/i{yi) is 
<5^<5(1 — yi) or the PDF of the parton j in the resolved photon i — 7 , and a(jk -x H + X) 
is the partonic cross section. I 11 NRQCD through relative order v 2 . the latter is factorized 
as 


a{ij^H + X) = Y J 


Fij(n ) 


h(0 H (n)) + 


Gij(n) 

m c ~nri c 


a O(n) 


h(F H (n)) 


a V(n) 


( 2 ) 


where O n ( n ) is the four-quark operator pertaining to the transition n —y H at leading order 
(LO) in v, with mass dimension do { n ); V H (n ) is related to its v 2 correction and carries mass 
dimension <ip( n ) = do( n ) + 2; and F lJ (n) and G t j(n) are the appropriate SDCs of the partonic 
subprocesses i + j —> cc(n ) + A". Working in the fixed-flavor-number scheme, the parton i 
runs over the gluon g and the light quarks q = u,d,s and antiquarks q. 

According to the velocity scaling rules 3], the leading contributions to direct J/-0 and 
ip' production are due to the 3 S'[^, 3 Sf\ ^q 8 ^, and 3 Pj 8 ^ channels, and those to direct XcJ 
production are due to the 3 P^ and 'iS'j 8 ' channels. Accordingly, prompt J/ijj photoproduction 
and hadroproduction proceeds at LO through the partonic subprocesses 

9 + 7 -> cc( 3 s!‘- 81 , 1 Si 81 , 3 Pj 8| )+g, 

9 ( 9)+7 -+ cc( 3 Sl s| , 1 4 8| , 3 Pi 81 ) + 9@. 

9 + 9 -+ ce( 3 si 1 'V4t 3 b 1,8| )+9, 

9®+9 -+ cc( 3 Sf 1 , 1 Sj, 8| , 8 Pi 1 ’ 81 ) + 9(9), 

9 + 9 -» cc( 3 Si S| , I 4 8 Pd 1 ' S ')+9- (3) 


We adopt the definitions of the relevant four-quark operators O h (n) from the literature 
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23] and define the corresponding four-quark operators V H (ri) as: 
V J/ ^( 3 Sf 1 ) = x + H.c., 


VWfS®) = x ] T a ^ JH ,aj/^T a 


- l -D) x + H.c., 


vJ/'/’C P [ f ] ) = x tcr * 


T^a^aj/^cr'T* (--0 


pXcO^pI 11 ) = axc0 )^t 


—U-O- 


-^) X + H.c., 


X + H.c., 


pxoi( 3 pW) = jLt /'-IS’ x <r) a Xcl )^ x cr) X + H.c., 


pxc 2 ( 3 pW) = * t /"_1S ( V J)N ) ^(4c 2 a Xc 2 )S (-7$) X + H.c., 


^( 3 4 81 ) = xVr>(4,a^VW-^Sj X + H.c., 

P^( 3 M 81 , 3 M 81 ) = y|x t a i T B ^(aW)^ t ^^ y r B X + H.c., 


(4) 


where = (SV J + S- 7 'cr i )/2 - S ■ cr^/3 and = (-i/2) 2 (‘S‘S' - S 2 ^/^ 
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Our definition of the S-D mixing operator 7 :>// ( 3 5'| 8 ^ 3 Z)^) differs from that in Refs, 
where a linear combination of ’P' H ( 3 S'f^) and P i? ( 3 <S'{ 8 ^ 3 -Df^) i n Eq. (J4]) is used instead. At 
order v 2 , there are no heavy-quark spin symmetries among the 0 H (n ) operators, but they 
still hold among the V H (n) operators, yielding the relationships 


(P J/ ^( spj 81 )) = 


{v^e sf 1 )) = 2J+1 
(v^epi 1 )} = 


2J + 1 
1 


(V^CP?)), 

(pxc( 3 5 [ 8l )) ; 

(pxcj ( 3 pj 1 ]) ) , 


(pxco^plsUpj 81 )) = -i(VXd^ s m3 D f ] )) = 2(V Xc2 Csf\ 3 Df ] )). (5) 

o 

The SDCs F tJ (n) and Gij{n ) may be obtained perturbatively by matching the QCD and 
NRQCD calculations via the condition 

^(cc)| p ert QCD = ^ d ”| n) 1 4 (° I G^(K) |0) | pert NRQCD + ZL-4 (O^n^A) 1°) I pert NRQCD- 


m c 


The left-hand side of 
developed in Ref. 


( 6 ) 

Iq. dnj) may be computed directly using the spinor projection method 


24J, by which the product of Dirac spinors v(P/2 — q)u(P/2 + q) is 
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projected onto the considered 2S+1 Lj state in a Lorentz-covariant form. In an arbitrary 
reference frame, the four-momenta P/2 + q and P/2 — q of the heavy quark and antiquark 
may be related to those in the quarkonium rest frame as 


P 


q = L 


Pr 


+ q , 


p 


— q — L 


<7 , 


(7) 


where P' r J = (2E q ,0), E q = \/m 2 c + q 2 , 2 q is the relative three-momentum between the 
two quarks in the quarkonium rest frame, and IX v is the Lorentz transformation matrix for 
the boost from the quarkonium rest frame to the considered reference frame. To all orders 
in v 2 , the projectors onto the spin-singlet (S = 0) and spin-triplet (S = 1) states in the 


quarkonium rest frame read 


25] 


E / . _ / . / 1 X | \ Eq T ‘TTlr 

v(—q, A 2 )w(<7, Ai) Ai; - , A 2 | 0 , Oy = —-^=— 


cl ■ q 

E 0 + m c 


x 7' 


,1 + 7 ° 


1 + 


Y v{-Q, A 2 )n(q, Ai 

Al,A2 


1 , 7 \ ! i \ _ Eg + m c 

, Ai; -, A 2 11, e ) — — —j =— I 1 — 


cl ■ q 

E q + m c 

cl ■ q 


7 


E 0 + m c 


x cl ■ e 


1 + 7° 


1 + 


In an arbitrary reference frame, they become 

Y A 2)h(g, Ai) /i, Ai; A 2 |0, o\ 


2^2 (Eg + m c ) V 2 


ol ■ q 

Eg + m c 


— Q-m. 
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( 8 ) 


x Y 


2 E n 


2 E n 


E »(-*.*>)«(«. AO (5. a.; 5, A,n, e ) = mE } +mc 


+ <j + m. 


m, 


x /- 


2 Eq 


2 Eq 


— + f + m c 


(9) 


Note that the normalization of the Dirac spinors is uu = — vv = m 2 . With the help of the 
spinor projection method, the partonic scattering amplitude M(ij —> cc(n) + X ) may then 
be expanded in the relative momentum q. To this end, we write 


TYl 

M(ij -> cc(n) + X) = -±A(q), 

V 


( 10 ) 
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where the factor a/ mn c /E q stems from the relativistic normalization of the cc(n) state and 


Ma) = E ^\S,sA{3,k;3,l\l(S,a)) 


Ai,A 2 k,l 

x A\ij -> c Al ,fc 


P 


+ q c X2 , 


p 


-q)+X 


( 11 ) 


ffere, (3, k] 3, Z|l) = 8ki/\f7P c an d (3, k; 3, l\8, a) = \/2 TJ? Z are the color-SU(3) Clebsch- 
Gordan coefficients for the cc(n) pair projected onto CS and CO states, respectively, and 
A(ij —> c Xl , k {P/2 + q)c\ 2 j(P /2 — q) + A") is the standard Feynman amplitude. Defining 

d N A(q) 


A-ai ■■■aN (d) 


dq ai ■ ■ ■ dq aN 
we may write the expansion of Eq. (TTT1) in q as 


( 12 ) 


q =0 


A{q) — A(0) + q ai A ai (0) + -q ai q a2 A aia2 ( 0) + -q ai q° 2 q a3 A aiol2a3 (0) + • • • . (13) 

2 o 

For S- and D-wave states, only the terms with even powers in q contribute, while for P-wave 
states it is the other way around. To calculate the relativistic corrections for the production 
of S- and P-wave states, we need to decompose the higher-rank tensor products of q factors 
in Eq. (TT3l) into their irreducible representations and to retain the L = S and L = P terms, 
respectively. Through order v 2 , we thus obtain, for n = 3 Sf\ 


rn r 


M(ij —» cc(n) + A") = . / —— 


E n 


A(Q) + \ALn a ^A aia2 (0) 
6 


(14) 


and, for n = 3 Pj^, 3 Pf\ 


M (ij —> cc(n ) + X) = 



C fti 

q 1 


4J0)- 


30 


TT Q2Q!3Q!4 A 

■ rL a20L3QL4 


( 0 ) 


(15) 


where IF 1 " 2 = -g a 102 + P ai P<* 2 / (4E 2 ) and IF 1 " 2 " 3 " 4 = IP 1 ® 2 IF 3 " 4 + IP 1 " 3 IF 2 " 4 + 
n“ 2 “ 3 n"i Q 4. In the D -wave case, we only need the amplitude at LO in v 2 , which is 


M(ij cc( 3 Df) + X) = ^q a2 A aia M. 

Il/q Z 


FEE 1 


(16) 


We are now in a position to perform the matching between the calculations in NRQCD 
and full QCD. We thus obtain 

Fij (n) 


m, 


l O(n) 


-4 


1 

2 s 


dLIPS|Mij(n)r, 


EM. = 1 / rfLIPS (K\Mij(n)\ 2 + |A r y(n)| 2 ), 


m. 


2s 


(V) 
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where yfs is the invariant mass of the incoming partons and 


| M ,»| 2 = J >( 0 )| : 

L z 

i%(«)i 2 = £ 1 0 

L, 


dq 2 


q 2 =0 

m c 

~En 


|A(0)| : 


+ -n QlQ2 Re[^(0)A Qia2 (0)] 


(18) 


q 2 =0 


forn = 3 5 , i 1] , 3 5{ 8l , 1 5S 8] ; 


|M,,(n)| 2 = ]T |(1, L z - 1,S Z \ J,J z )e*M a (0)f 

S Z ,L Z ,J Z 

mnt h rm " 


q 2 =0 


* S Z ,W; 

1 


E„ 


|(1, L z - 1,^| J,J z )e* L a A a (0)\' 


15 


Y (1, L z ; 1, S Z \J, J z ) Y (1,^1; l,S' g \J, J' z ) 


S Z ,L Z ,J Z 


^ j-j aiQ2Q;3 




(19) 


q 2 = 0 


for n = 3 Pj^ , 3 Pf ] ; and 


\Nij(n)\ 2 = Y (^L z -l } S z \l,J z )Re[eTfA*(0)A a ^0)] 

S Z ,L Z ,J Z 


( 20 ) 


q 2 =0 


for the 3 S'^- 3 Zlf^ mixing term. Here, e^ z (e^f) is the polarization four-vector (four-tensor) 
for L = P (D), the symbol 22 also implies the summation over the polarizations of the other 
external partons, and the bar implies the average over the spins and colors of the incoming 
partons and those of the ce(ri) state. The factor K in Eq. (fI7|) contains the v 2 corrections 
to the phase space and depends on the kinematic variables that we are interested in. In 
the cases under consideration here, we have K = —4/(s — 4m 2 ). We generate the Feynman 


diagrams using the FeynArts package 26] and compute the amplitude squares using the 


FeynCalc package [27]. In the Appendix, we present our results for \Nij(n)\ 2 . 

We reproduce the well-known results for |Mjj(n)| 2 , which were called |^4| 2 in Ref. 28]] and 
\M'\ 2 in Ref. 29]. We may also compare some of our results for \Nij{n)\ 2 with the literature. 
The v 2 corrections to direct J/i[} photoproduction were first studied in Ref. 30] within the 
relativistic quark model, which accounts for the CS contributions. We reproduce Eq. (23) 


in Ref. 30] if we do not expand E q , but set E q = mj/,/,/2. The relativistic corrections to 
direct J/ijj hadroproduction were considered in Ref. 19j within NRQCD. We find agreement 
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(a) (b) (c) 


FIG. 1: Ratios R(n) for J/'ij) direct photoproduction under HERA II kinematic conditions as 
functions of (a) p\, (b) W, and (c) 2 . 


with Ref. 19j, except for some typographical errors in Eq. (A4) therein, which corresponds 
to |A^( 3 Sf ] )| 2 in our notation. Furthermore, the 3 S'[ 8 ^- 3 d[ 8 ^ mixing contribution was not 
considered there. Apart from correcting the misprints in Eq. (A4) of Ref. 19|, our paper 
reaches beyond the available literature by studying the CO contributions to direct J/^j 
photoproduction, the 3 A| 8 - 3 d| 8 ' mixing contribution to direct J/-0 hadroproduction, and the 
feed-down contributions to J/'ip hadroproduction. 


III. PHENOMENOLOGICAL RESULTS 


We are now in a position to investigate the phenomenological significance of the v 2 cor¬ 
rections in prompt J/ip photoproduction and hadroproduction. In our numerical analysis, 
we use m c = 1.5 GeV, a = 1/137.036, the LO formula for a " 1 (ji r ) with rif = 4 active quark 
flavors and asymptotic scale parameter Aq^, d = 215 MeV 31], the CTEQ6L1 set for proton 
^DFs 31], the photon flux function given in Eq. (5) of Ref. 32] with Q^ax = 2.5 GeV 2 
33], and the choice p r = /i/ = + 4m 2 for the renormalization and factorization scales. 

According to Eqs. (JT|) and (J2J), the hadronic cross sections of the direct photoproduction and 
hadroproduction of the charmonia H = J/ijj , Xcj , ^, differential in some observable x , may 
be generically written as 


da 

dx 


E 


dF(n) (O h (7i)} dG(n) ( V u {n)) 


dx 


771 . 


do(n)—^ 


+ 


dx 


m, 


dv(n)~^ 


( 21 ) 


where it is understood that dF{n)/dx = 0 if n stands for 3 S'j 8 '- 3 Zl| 8 ^ mixing. For 
all other channels n, the relative v 2 corrections are R(n)(V H (n))/(ml(0 H (n))), where 
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R(n) = dG(n)/dx/(dF(n)/dx ) is a dimensionless ratio. To standardize the numerical 
discussion, we also introduce the quantity R(n) = dG{n)/dx/(dF( 3 sf^)/dx) for the case 
when n stands for 3 5'[ 8 - 3 Dj 8 ^ mixing. According to the velocity scaling rules [3], the weight 
(V H (n))/(m 2 (Q H (n))) of R(n) is of order v 2 , which is approximately 0.23 for charmonium 
34,135j]. For definiteness, we ignore these weights in the following and concentrate on the 
SDC ratios R(n) instead. Furthermore, we limit ourselves to the direct production of J/ip 
mesons and their production via the feed down from the \ c j mesons. In the latter case, 
the branching fractions Bjcj —>■ J /V* + -V) drop out in the ratios R(n). However, there 
remains a kinematic effect on the transverse momentum. Since pt = p'r^ M XcJ — Mj/^, 
we may approximate p^ v ’ = p^ cj Mj/^/ M Xc „ with Mj/$ = 3.097 GeV, M Xc0 = 3.415 GeV, 
M Xcl = 3.511 GeV, and M Xc2 = 3.556 GeV 36]. 


We consider three typical experimental environments, namely, Run II at HERA, Run I at 
the Tevatron, and the LHCb setup at the LHC. At HERA II, the cross section of prompt J/%) 
photoproduction was measured at center-of-mass energy \fS = 319 GeV differential in p?, 
W = a/( p p + p 7 ) 2 , and 0 = pj/p-pp/Pj-Pp 33], where p p , p 7 , and pj are the four-momenta 
of the proton, photon, and J/ip meson, respectively, imposing in turn two of the acceptance 
cuts 1 GeV 2 < 60 GeV < W < 240 GeV, and 0.3 < z < 0.9. In Figs. [Ela)-(c), the 

Py, W, and z distributions of R(n) are shown for n = :i S \ 1 ', 3 5'j 8 ^, , 3 Pj 8 ', and 3 S'[ S '- 3 i4^ 

mixing. We recall that n = 3 P^ is prohibited at this order. For simplicity, we ignore resolved 
photoproduction, the SDCs of which also contribute to hadroproduction to be studied below. 
In Run I at the Tevatron, the pr distribution of prompt J Jib hadroproduction was measured 
at \/~S = 1.8 TeV in the pseudorapidity range |p| < 0.6 37], and it was measured by LHCb 
at \fS = 7 TeV in the rapidity range 2.0 < y < 4.5 38]. In the latter two cases, we 
exclude the region pr < 3 GeV, where the application of fixed-order perturbation theory 
is problematic. In Figs. [2] and [3j the Pt distributions of R(n) are shown for the Tevatron 
and the LHC, respectively. In each figure, part (a) refers to n = 3 5'[^, 3 5 , [ 8 ^, , and 

3 Sf ] - 3 Df ] mixing in direct J/if) production, and part (b) refers to n — :i Pj in the feed down 
from the respective \ c j meson. In part (b), we do not show R(n) for n = 3 S'[ 8 ^ and 3 S'[ S ^- 3 Z1^ 
mixing because these results emerge from their counterparts in part (a) just by rescaling the 
Pt axis as explained above. In a similar way, part (a) carries over to the feed down from the 
J meson. 

We observe from Figs. [£H2]that the ratios R(n) may be throughout positive or negative 
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(a) (b) 


FIG. 2: Ratios R(n) for J/^ hadroproduction (a) in the direct mode and (b) via the feed down 
from XcJ mesons under Tevatron Run I kinematic conditions as functions of px- 


or change sign as x is varied and typically have a magnitude of order unity or larger. In fact, 
-RpSj 8 ]) for hadroproduction reaches the value —3 at px = 3 GeV, as may be seen in Figs. [2] 
and[3j An exception to this rule is /?,(' 3 5'[^) for hadroproduction, which is roughly one order of 
magnitude smaller. Because of the above estimations of the weights (V H (n)) / (m 2 ((D H (n))) 
and v 2 , R(n ) = 2 is likely to imply a v 2 correction of some 50% in the channel n. However, it 
is unlikely that the inclusion of the v 2 corrections would turn a given hadronic cross section 
negative, the more so as the v 2 corrections in the various channels n may be of either sign thus 
allowing for cancellations. Comparing Figs. HHSl we observe that the relative importance 
and the px dependencies of the various ratios R{n ) greatly differ between photoproduction 
and hadroproduction. If a fit is constrained to data for which the contributing ratios R(n ) 
are approximately independent of x, as in the case of hadroproduction for px > 7 GeV say, 
then the inclusion of the v 2 corrections would not improve the quality of the fit, but rather 
induce strong correlations between ( 0 H (n )) and (' P H {n )) for each n, both as for sign and 
magnitude. To reduce these correlations or even enable independent determinations of the 
LDMEs in the enlarged set, it is indispensable to reduce the low-p^ cut in hadroproduction, 
to px > 3 GeV say, and to include photoproduction data in the fit. In a fit with such a low- 
px cut in hadroproduction, the inclusion of the v 2 corrections might enhance the goodness 
because the notorious lack of turnover of the NLO NRQCD predictions in the low-py region, 
which causes the latter to overshoot the experimental data there, might be cured at least 
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(a) (b) 


FIG. 3: Ratios R(n) for J/^ hadroproduction (a) in the direct mode and (b) via the feed down 
from XcJ mesons under LHCb kinematic conditions as functions of px- 

to some extent by the fact that all the ratios R{n ) are amplified as px is decreased. In 
conclusion, the addition of v 2 corrections is expected to have a significant impact on state- 
of-the-art determinations of CO LDMEs through global fits to prompt J/^ production data 
and might improve the goodness of such fits. 


IV. CONCLUSIONS 

We performed a systematic study of the v 2 corrections to the yields of prompt J/ijj pho¬ 
toproduction and hadroproduction, providing the relevant SDCs in analytic form. Specif¬ 
ically, this includes the partonic subprocesses listed in Eq. (|3]) in combination with the cc 
Fock states n = ?J S\ 1 ^, 3 Pj'^, 3 S'[ 8 ^, , 3 Pj 8 ^, and 3 S'[ 8 ^- 3 Zl[ 8 ' mixing. We compared our results 

with the literature as far as the latter goes. We assessed the phenomenological significance 
of the v 2 corrections in the various channels by studying their ratios with respect to the 
corresponding LO results. Assuming the relevant LDMEs (O h ( n)) and (' P H (n )) to obey 
the hierarchy predicted by the velocity scaling rules fj, we found that v 2 corrections of up 
to 50% are realistic. We thus conclude that it is indispensable to include v 2 corrections in 
determinations of CO LDMEs by global data fits, the more so as the v 2 corrections in the 
various channels greatly differ between photoproduction and hadroproduction. 
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V. APPENDIX 

In this appendix, we present analytic expressions for the nonvanishing SDCs \N i; j(n)\ 2 
appropriate for prompt J /-0 photoporduction and hadropro duct ion evaluated according to 
Eqs. (USD-®. For the partonic subprocess i(pi) + jijp-i) —* cc(P) + X, the Mandelstam 
variables are defined as s = (pi + p 2 ) 2 , t = {p\ — P) 2 \ q= o, and u = (p 2 — P) 2 \ q =o and satisfy 
s + t + u = 4m 2 . Note that our results for hadroproduction can also be applied to resolved 
photoproduction. 

A. Photoproduction 

9 + 7 -> ccpsl 11 ) +g: 


|Jv„( 3 si‘')l 2 = 


81927r 3 aa^ 


729 m c (s — 4 m;?) 3 (4m;? — t) 3 (s + t) 3 


x 


2048m 30 (3s 2 + 2 st + 3 1 2 ) - 256m® (5s 3 - 8 s 2 t + 4sf 2 + 5 1 3 ) - 64m°(15s 4 + 68 s 3 t 
+62 s 2 t 2 + 32 st 3 + 15f 4 ) + 16m 4 (21s 5 + 87 s 4 t + 130s 3 t 2 + 106s 2 t 3 + 51st 4 + 21t 5 ) 
—4m 2 (7s 6 + 30s 5 f + 59 s A t 2 + 64s 3 f 3 + 47s 2 t 4 + 18sf 5 + 7f 6 ) — st(s + t) (s 2 + st + f 2 )' 


( 22 ) 


9 + 7 -A cc( 3 Ff ] ) + g: 


g + i -+ ccO-S? 1 ) + g- 


lAW^f 1 )! 2 = tIA'„( 3 s! 1| )I 


15, 


/3 ebb 12 


(23) 
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pv s7 ( 1 sf , )| 2 =- 3 12S *Y- st - 

9 ml (Ami ~ s ) 3 (4 ml — t ) 3 (s + t ) 3 (—Ami + s + t) 
x [l31072m 44 - 4096m 42 (23s + lit) + 1024m 40 (33s 2 + 34st - 3t 2 ) - 256m 8 (s + t) 

(7s 2 - list - 29t 2 ) - 64 m 6 (l3s 4 + 94s 3 t + 170s 2 t 2 + 142st 3 + 37t 4 ) + 16m 4 (23s 5 
+112s 4 t + 214s 3 t 2 + 226s 2 t 3 + 112st 4 + 231 5 ) - 4m 2 (9s 6 + 50s 5 t + 116s 4 t 2 + 138s¥ 
+104s 2 t 4 + 38st 5 + 9t 6 ) + 5st(s + t) (s 2 + st + t 2 ) 2 (‘ 

9 + 7 -> cc( 3 Pj 8] ) + g: 


/3p[8K.2 = _ I 28 n 3 aa 2 s _ 

91 J 15 ml (s — 4m ;?) 4 (t — 4m ;?) 4 (s + t ) 4 (—4m;? + s + t) 

x [3145728m 48 (s + t) (s 2 + 8 st + t 2 ) - 65536m 46 (18s 4 + 113s 3 t + 370s 2 t 2 - 27st 3 
+18t 4 ) - 32768m 44 (39s 5 + 318s 4 t + 459s 3 t 2 + 929s 2 t 3 + 548st 4 + 39t 5 ) + 16384m 42 
(63s 6 + 519s 5 t + 1200s 4 t 2 + 2043s 3 t 3 + 1805s 2 t 4 + 629st 5 + 63t 6 ) - 2048m 40 (144s 7 
+1231s 6 t + 3673s 5 t 2 + 7074s 4 t 3 + 7874s 3 t 4 + 4843s 2 t 5 + 1301st 6 + 144t 7 ) + 512 m 8 
(75s 8 + 662s 7 t + 2626s 6 t 2 + 5981s 5 t 3 + 7990s 4 t 4 + 6371s 3 t 5 + 3206s 2 t 6 + 642st 7 + 751 8 ) 
— 128m 6 (15s 9 + 162s 8 t + 964s 7 t 2 + 2725s 6 t 3 + 4364s 5 t 4 + 4034s 4 t 5 + 2635s 3 t 6 + 1114s 2 t 7 
+132st 8 + 15t 9 ) + 16m 4 st(27s 8 + 502s 7 t + 1968s 6 t 2 + 3774s 5 t 3 + 4338s 4 t 4 + 3114s 3 t 5 
+1668s 2 t 6 + 502st 7 + 271 8 ) - 16m 2 s 2 t 2 (s + t)(25s 6 + 118s 5 t + 259s 4 t 2 + 310s 3 t 3 
+244s 2 t 4 + 103st 5 + 25t 6 ) + 31s 3 t 3 (s + t ) 2 (s 2 + st + t 2 ) 2 1 (25) 

9 + 7 -t cc( 3 S [ i\ 3 Df ] ) + g: 


]N„esf\ 3 Df ] )] 2 = 


2048^5/3 7r 3 aa 2 


■ 32768m 44 (s + t)(3s 2 


^ “ 81m c (s - 4m ;?) 4 (t - Am 2 f (s + t ) 4 L ^ ^ 

+5st + 31 2 ) — 4096m 42 (s 2 + 5st + t 2 ) (lls 2 + 18st + lit 2 ) — 1024m 40 (s + t)(10s 4 

—109s 3 t - 202sV - 109st 3 + 10t 4 ) + 256m 8 (36s 6 - 27s 5 t - 395s 4 t 2 - 640s 3 t 3 

—395s 2 t 4 - 27st 5 + 36t 6 ) - 64m 6 (s + t)(28s 6 - 3s 5 t - 250s 4 t 2 - 400s 3 t 3 - 250s 2 t 4 

-3st 5 + 281 6 ) + 16m 4 (7s 8 + 4s 7 t - 100s 6 t 2 - 323s 5 t 3 - 440s 4 t 4 - 323s s t 5 - 100s¥ 


+4st' + 71 8 ) + 4m 2 st(s + t) ( 6 s 6 + 22s 5 t + 51s 4 t 2 + 66 s 3 t 3 + 51s 2 t 4 + 22st 5 + 6 t 6 ) 

—s 2 t 2 (s + t ) 2 (s 2 + st + t 2 )^ (26) 
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?(?) + 7 -> cc( 3 Si 81 ) + q(q): 


|A r 9 (5h( 3 Sf')l' 


87r 3 aa 2 e 2 


x 


27 m b c st (4 ml — s ) 

640m® — 160m^(2s + t) + 4m 2 (27s 2 + lOst + 5t 2 ) — 11s (s 2 + t 2 ) 


(27) 


q(q) + 7 -> ccf'slf 1 ) + 9ta) : 


|A r ,(,-h(‘5f l )| 2 


2567r 3 aa 2 [m 2 (44s 3 + 92s 2 t — 4st 2 + 44t 3 ) — 5s(s + t ) (s 2 + t 2 )] , 

81m 3 (s — 4m;?) (s + t) 3 (—4 m;? + s + t) 


q(q) + 7 -> ecfi 3 ] 81 ) + q(q): 


i^,«h( 3 7 > l i| )r= 


2567r 3 aa 2 


X 


135m 3 (s — 4 m;?) (s + t) 4 (—4m;? + s + t) 

512 ml (5s 2 + 26st + 25t 2 ) + 64m* (s 3 - 23 s 2 t - 111st 2 - 19t 3 ) + 4m 2 (s + t) 


(57s 3 + 169s 2 t - 3 st 2 + 611 3 ) - 31s(s + t ) 2 (s 2 + t 2 ) 
g(g) + 7 -)• cc{*S?M?) + ?(?): 


(29) 


|/V s( 5 h ( 3 sf , , 3 c [ 8| )| 2 


16A/5/37r 3 aa 2 e 2 [32m 4 — 8m 2 (s + t) + s 2 + t 2 
9 m b c st 


(30) 


B. Hadroproduction 


q + g —> cc( 3 5{ 1] ) + q: 


ia , »( 3 7 1| )i 2 


15a s 

128a 


|A r „( 3 5|‘ l )| 2 


(31) 


g + g -> of/f 1 ) + S' 
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I_/v ( 3 pM')| 2 —_ fts _ 

Abm^st (4m% — s) 5 (4ml — t) 5 (s + t) 5 (—Am^ + s + t) 
x 251658240(s + t) 3 (lls 2 + 30ts + lit 2 ) m 24 - 4194304(s + t) 2 (5s + 3 1) (231s 3 + 847ts 2 

+921t 2 s + 385t 3 )m 22 + 1048576(s + t) 2 (4155s 5 + 19680ts 4 + 36066t 2 s 3 + 34246tV 
+17160t 4 s + 3795t 5 )m 20 - 262144(9735s 8 + 69804ts 7 + 219507t 2 s 6 + 399776t 3 s 5 
+466128t 4 s 4 + 359956t s s 3 + 181147t 6 s 2 + 54864t 7 s + 7755t 8 )m 18 + 65536 (15750s 9 
+119943ts 8 + 411033tV + 838246t 3 s 6 + 1131460t 4 s 5 + 1057300t 5 s 4 + 688706t e s 3 
+304133t 7 s 2 + 83403t 8 s + 10890t 9 )m 16 - 32768 (8955s 10 + 72307ts 9 + 267272t 2 s 8 
+600188t 3 s 7 + 913993t 4 s 6 + 992452tV + 782713t 6 s 4 + 445538t 7 s 3 + 177082tV 
+44767t 9 s + 5445t 10 )m 14 + 4096(l4235s n + 121663ts 10 + 482477t 2 s 9 + 1181531t 3 s 8 
+1997582t 4 s 7 + 2459484t 5 s 6 + 2258884t 6 s 5 + 1553862t 7 s 4 + 791411t 8 s 3 + 287757t 9 s 2 
+67783t 10 s + 7755t n )m 12 - 1024(7575s 12 + 69080ts n + 295378t 2 s 10 + 788288t 3 s 9 
+1469633t 4 s 8 + 2024628t 5 s 7 + 2117532t 6 s 6 + 1697668t 7 s 5 + 1040113t 8 s 4 + 480268t 9 s 3 
+160638t 10 s 2 + 35180t n s + 3795t 12 )m 10 + 256 (2415s 13 + 24464ts 12 + 115885t 2 s n 
+342339tV° + 708615t 4 s 9 + 1090700t 5 s 8 + 1287226t 6 s 7 + 1180426t 7 s 6 + 842220tV 
+464375t 9 s 4 + 194239t 10 s 3 + 58985t n s 2 + 11684t 12 s + 1155t 13 )m 8 - 64 (345s 14 
+4528ts 13 + 25895t 2 s 12 + 89188t 3 s n + 211442t 4 s 10 + 369630t 5 s 9 + 495006t 6 s 8 + 517164t 7 s 7 
+424126t 8 s 6 + 271810t 9 s 5 + 134322t 10 s 4 + 49688t n s 3 + 12995t 12 s 2 + 2128t 13 s + 165t 14 )m° 
+16st(s + t) (259s 12 + 2323ts n + 9736t 2 s 10 + 25649t 3 s 9 + 47816t 4 s 8 + 66583t 5 s 7 + 71232t 6 s 6 
+58883t 7 s 5 + 37456t 8 s 4 + 17829t 9 s 3 + 6056t 10 s 2 + 1303t n s + 139t 12 )m 4 - 4s 2 t 2 (s + t ) 2 
(s 2 + ts + t 2 ) 2 (87s 6 + 372ts 5 + 746t 2 s 4 + 854t 3 s 3 + 666 t 4 s 2 + 292t 5 s + 67t 6 ) m\ 

+13s 3 t 3 (s + t) 3 (s 2 + ts + t 2 ) 4 (32) 

9 + g ->• cc( 3 Pj 1] ) + g: 


16 



|A^( 3 Pi [1] )| 2 =----- 

45m® (s — 4mj?) (4m^ —t) 5 (s + t ) 5 

x 2097152m 20 (s + tf (s 2 + t 2 ) - 131072m 18 (lls 5 - 35s 4 t - 40sV + 5st 4 + lit 5 ) 
-32768m 46 (125s 6 + 864s 5 t + 1419s 4 t 2 + 1164s 3 t 3 + 739s 2 t 4 + 404st 5 + 851 6 ) 
+8192m 14 (818s 7 + 4579s 6 t + 9054s 5 t 2 + 9591s 4 t 3 + 7291s 3 t 4 + 4914s 2 t 5 + 2379st 6 
+538t 7 ) - 2048m 12 (2106s 8 + 11940s 7 t + 27483s 6 t 2 + 35428s'¥ + 31746s 4 t 4 + 23628 s 3 t 5 
+14863s 2 t 6 + 6300st 7 + 1306t 8 ) + 512m 10 (2905s 9 + 17675s 8 t + 46718s 7 t 2 + 71799s 6 t 3 
+75989s 5 t 4 + 63429s 4 t 5 + 44939s 3 t 6 + 25058s 2 t 7 + 9275st 8 + 1705t 9 ) 

— 128 m 8 (2265s 10 + 15280s 9 t + 46309s 8 t 2 + 84008sV + 105380s 6 t 4 + 101476s 5 t 5 
+79720s 4 t 6 + 51188s 3 t 7 + 24729s 2 * 8 + 7940st 9 + 1265t 10 ) + 32m 6 (944s n + 7381s 10 t 
+26060s 9 * 2 + 56013s 8 t 3 + 83898s 7 * 4 + 95046s 6 * 5 + 85066s 5 * 6 + 61038s 4 * 7 + 34133s 3 * 8 
+14020s 2 * 9 + 3821st 10 + 504* 11 ) - 8m 4 (164s 12 + 1700s 11 * + 7455s 10 * 2 + 19660s 9 * 3 
+35909s 8 * 4 + 48760s 7 * 5 + 51084s 6 * 6 + 41660s 5 * 7 + 26329s 4 * 8 + 12440s 3 * 9 + 4175s 2 * 10 
+900s* n + 84t 12 ) + 2m 2 s*(s + *) (s 2 + st + t 2 ) 2 (106s 6 + 490s 5 * + 981s 4 * 2 + 1092s 3 * 3 
+821s 2 * 4 + 330s* 5 + 66 t 6 ) - lls 2 t 2 (s + t ) 2 (s 2 + st + t 2 ) 4 (33) 

g + g -»• cc( 3 Pj 1] ) + g: 
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\N gg ( 3 P^ ] )\ 2 = --- 5 - 10^4 - 

225m® st (Ami ~ s ) 5 (4ml — t)° (s + t) 5 (—Ami + s + t) 

x 100663296(s + t) 3 (l5s 2 + 38 ts + 15t 2 ) m 24 - 8388608(s + t) 2 (315s 4 + 1312ts 3 
+2010tV + 1252t 3 s + 315t 4 )m 22 + 524288(s + t) (4620s 6 + 25287ts 5 + 58998t 2 s 4 
+73262tV + 53438tV + 21927t 5 s + 4140t 6 )m 2 ° - 262144(5550s 8 + 36807ts 7 
+109272tV + 189017t 3 s 5 + 212562t 4 s 4 + 162877t 5 s 3 + 84442tV + 27387 t 7 s 
+4230t 8 )m 48 + 32768(18360s 9 + 126189ts 8 + 389652tV + 714467t 3 s 6 + 879782tV 
+782522tV + 517507t 6 s 3 + 248752tV + 78429t 8 s + 11880t 9 )m 46 - 8192(21240s 10 
+155344ts 9 + 504167tV + 967724tV + 1252087t 4 s 6 + 1199524tV + 911047tV 
+558524tV + 262447t 8 s 2 + 81544t 9 s + 11880t 10 )m 44 + 2048(17100s n + 136829ts 10 
+477883tV + 972667 t 3 s 8 + 1311058tV + 1293789t 5 s 6 + 1038389t 6 s 5 + 732018tV 
+446587t s s 3 + 210243tV + 63149t 10 s + 8460t n )m 42 - 512(9180s 12 + 82966ts n 
+321992t 2 s 10 + 716914t 3 s 9 + 1033387t 4 s 8 + 1055592tV + 855762 t 6 s 6 + 634352tV 
+454727tV + 283274tV + 129252t 1( V + 35566t n s + 4140t 12 )m 4 ° + 128(2940s 13 
+31855ts 12 + 143846t 2 s n + 365607tV° + 592170tV + 663334t 5 s 8 + 565754tV 
+427154t 7 s 6 + 326754t s s 5 + 239410tV + 141047t 10 s 3 + 57346t n s 2 + 13375t 12 s 
+1260t 13 )m 8 - 32(420s 14 + 6548ts 13 + 37651tV 2 + 116096t 3 s n + 224575tV° 
+299040tV + 299220 t 6 s 8 + 251484tV + 199880tV + 152480tV + 99555t 10 s 4 
+48976t n s 3 + 16051t 12 s 2 + 2828t 13 s + 180t 14 )m 6 + 8 st(s + t)(464s 12 + 4070ts n 
+15035tV° + 32905t 3 s 9 + 49483tV + 56786 t 5 s 7 + 54168t 6 s 6 + 44806t 7 s 5 + 32183t 8 s 4 
+18405tV + 7675t 10 s 2 + 2030t n s + 224t 12 )m 4 - 2 s 2 t 2 (s + t ) 2 (s 2 + ts + t 2 ) 2 (138s 6 
+552ts 5 + 943tV + 922t 3 s 3 + 783t 4 s 2 + 392t 5 s + 98t 6 )m 2 + 7 s 3 t 3 (s + t) 3 (s 2 + ts + t 2 ) 4 

(34) 


g + g ->• cc( 3 Sf ] ) + g: 
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|A^ 9 ( 3 si 8l )| 2 =---^- 3 - 

54m 3 (s — 4 m%) (4 -m% — t ) 3 (s + t ) 3 

x 16384m 34 (87s 2 + 22 st + 87 1 2 ) - 4096m 32 (l4s 3 - 449s 2 t - 221st 2 + 14t 3 ) 

-2048m 30 (480s 4 + 1969s 3 t + 2182s 2 t 2 + 1384st 3 + 399t 4 ) + 256m 8 (2910s 5 
+11616s 4 t + 17509s 3 t 2 + 15433s 2 t 3 + 8178st 4 + 2100t 5 ) - 64m 6 (4048s 6 + 17148s 5 t 
+31739s 4 t 2 + 35722s 3 t 3 + 25841s 2 t 4 + 11412st 5 + 2590t 6 ) + 16m 4 (2754s 7 + 12968s 6 t 
+28779s 5 t 2 + 39352s 4 t 3 + 35950s 3 t 4 + 22137s 2 t 5 + 8270st 6 + 1620t 7 ) 

-108m 2 (s 2 + st + t 2 ) (27s 6 + 131s 5 t + 269sV + 302s 3 t 3 + 227s 2 t 4 + 89 st 5 + 15t 6 ) 
+297st(s + t) (s 2 + st + t 2 ) 3 (35) 

g + g -»■ cc( 3 <S(j 8] ) + g~- 

|7v gg ( 3 4 81 )| 2 = - 3 -- 

3m 3 st (4m;? — s ) 3 (4m;? — t ) 3 (s + t) 3 (—4m;? + s + t) 
x 65536m 36 (5s 3 + 20s 2 t + 8 st 2 + 51 3 ) - 16384m 34 (25s 4 + 116s 3 t + 110s 2 t 2 + 56st 3 

+25t 4 ) + 12288m 32 (28s 5 + 129s 4 t + 189s 3 t 2 + 137s 2 t 3 + 65st 4 + 20t 5 ) 

-2048m 30 (87s 6 + 426s 5 t + 801s 4 t 2 + 806s 3 t 3 + 501s 2 * 4 + 204st 5 + 45t 6 ) 

+256m 8 (222s 7 + 1182s 6 t + 2653s 5 t 2 + 3385s 4 t 3 + 2749s 3 t 4 + 1501s 2 t 5 + 522st 6 + 90t 7 ) 
—64 m 6 (180 s 8 + 1020s 7 t + 2591s 6 t 2 + 3964s 5 t 3 + 4006s 4 t 4 + 2812s 3 t 5 + 1355s 2 t 6 
+408st 7 + 60t 8 ) + 16 m 4 (85s 9 + 516s 8 t + 1456sV + 2561s 6 t 3 + 3096s'V + 2688s 4 t 5 
+1673s 3 t 6 + 724sV + 192st 8 + 251 9 ) - 4m 2 (s 2 + st + t 2 ) 2 (17s 6 + 88 sH + 167 sH 2 
+168s 3 t 3 + 119s 2 t 4 + 40st 5 + 51 6 ) + 5st(s + t) (s 2 + st + t 2 ) 4 (36) 

g + g^ cc( 3 Pf ] ) + g: 
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IN f 3 pl 8 h| 2 =_ 71 Qs _ 

m b c st (s — Ami) 41 (t — Ami) 4 ( s + ¥ ( — 4 ml + s + t) 

x 1048576m 20 (ll5s 4 + 435s 3 ! + 564s 2 ! 2 + 295s! 3 + 115! 4 ) - 524288m 48 (345s 5 
+1571s 4 ! + 2683s 3 ! 2 + 2148s 2 ! 3 + 1156s! 4 + 345! 5 ) + 65536m 16 (2235s 6 + 11102s 5 ! 
+22740s 4 ! 2 + 23720s 3 ! 3 + 16400s 2 ! 4 + 7932s! 5 + 1955! 6 ) - 16384m 14 (4710s 7 
+24932s 6 ! + 57177s 5 ! 2 + 71731s 4 ! 3 + 60491s 3 ! 4 + 37287s 2 ! 5 + 16102s! 6 + 3450! 7 ) 
+4096m 12 (6660s 8 + 38018s 7 ! + 94927s 6 ! 2 + 134978s 5 ! 3 + 132094s 4 ! 4 + 97708s 3 ! 5 
+54647s 2 ! 6 + 21708s! 7 + 4140! 8 ) - 1024m 10 (6390s 9 + 39787s 8 ! + 107560s 7 ! 2 
+168509s 6 ! 3 + 183028s 5 ! 4 + 154598s 4 ! 5 + 104319s 3 ! 6 + 54540s 2 ! 7 + 20307s! 8 + 3450! 9 ) 
+256m 8 (4055s 10 + 27855s 9 ! + 82334s 8 ! 2 + 142320s 7 ! 3 + 169197s 6 ! 4 + 156936s 5 ! 5 
+120977s 4 ! 6 + 77050s 3 ! 7 + 38034s 2 ! 8 + 13095s! 9 + 1955! 10 ) - 64m 6 (1530s n 
+12009s 10 ! + 40145s 9 ! 2 + 79031s 8 ! 3 + 106482s 7 ! 4 + 110003s 6 ! 5 + 93903s 5 ! 6 + 68322s 4 ! 7 
+40531s 3 ! 8 + 18105s 2 ! 9 + 5449s ! 10 + 690! 11 ) + 16m 4 (255s 12 + 2610s 11 ! + 10798s 10 ! 2 
+26205s 9 ! 3 + 43547s 8 ! 4 + 54320s 7 ! 5 + 53906s 6 ! 6 + 43520s 5 ! 7 + 28607s 4 ! 8 + 14485s 3 ! 9 
+5278s 2 ! 10 + 1230s! 11 + 115! 12 ) - 4m 2 s!(s + !) (s 2 + s! + ! 2 ) 2 (150s 6 + 696s 5 ! 

+1299s 4 ! 2 + 1358s 3 ! 3 + 1059s 2 ! 4 + 456s! 5 + 90! 6 ) + 31s 2 ! 2 (s + !) 2 (s 2 + s! + ! 2 ) 4 1 (37) 

9 + 9 ->• cc( 3 5j 8] , 3 D [ * ] ) + g: 
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| N gg ( 3 4 81 , 3 Df ] )\ 2 = -=-- 4 - 

108\/l5 ml (s — 4 ml) 4 (* — 4 m%) (s + *) 4 

x [ - 524288m 20 (4463s 2 + 25840s* + 5057* 2 ) + 65536m 38 (s + *)(7767s 2 + 145760s* 
+13689* 2 ) + 16384m 36 (40922s 4 + 83203s 3 * - 78401s 2 * 2 + 46753s* 3 + 25451* 4 ) 

-4096m 34 (s + *) (41543s 4 + 579922s 3 * + 620206s 2 * 2 + 525220s* 3 + 20753* 4 ) - 3072m 32 
(40919s 6 - 153453s 5 * - 767319s 4 * 2 - 1075115s 3 * 3 - 737820s 2 * 4 - 133812s* 5 + 45008* 6 ) 
+256m 30 (s + *)(302824s 6 + 248112s 5 * - 1111219sV - 1766115s 3 * 3 - 1102093s 2 * 4 
+250353s* 5 + 304390* 6 ) - 64m 8 (287151s 8 + 1147982s 7 * + 1476785s 6 * 2 + 462324s 5 * 3 
-288201s 4 * 4 + 391881s 3 * 5 + 1418006s 2 * 6 + 1132457s* 7 + 286503* 8 ) + 16m 6 (s + *) 
(136080s 8 + 645678s 7 * + 1174264s 6 * 2 + 1234125s 5 * 3 + 1166106s 4 * 4 + 1190394s 3 * 5 
+1143268s 2 * 6 + 640305s* 7 + 136080* 8 ) - 4m 4 (27216s 30 + 226800s 9 * + 774349s 8 * 2 
+1531962s 7 * 3 + 2158036s 6 * 4 + 2377567s 5 * 5 + 2135518s 4 * 6 + 1512603s 3 * 7 + 768949s 2 * 8 
+226800s* 9 + 27216* 30 ) + m 2 s*(s + *) (18144s 8 + 104976s 7 * + 274390s 6 * 2 + 449444s 5 * 3 
+528929s 4 * 4 + 449390s 3 * 5 + 274363s 2 * 6 + 104976s* 7 + 18144* 8 ) 

—1620s 2 * 2 (s + *) 2 (s 2 + s* + * 2 ) 3 (38) 

q(q) + g ^ cc( 3 P^ 1] ) + q(q): 


\ N g(q)gC p o ] )\ 2 = -- 5 (12m 2 - *) [87040m 30 - 256m 8 

' q(q)9K ° 405m c 5 * (4m2 - s) (4m2 - *) 5 V c 

(285s + 278*) + 256m 6 (115s 2 + 166s* + 81* 2 ) - 32m 4 (ll5s 3 + 311s 2 * + 243s* 2 + 77* 3 ) 
+4m 2 * (96s 3 + 188s 2 * + 120s * 2 + 23* 3 ) - 13s* 2 (2s 2 + 2s* + * 2 ) 1 (39 

q(q) + g -»• cc( 3 pj 1] ) + g(q): 


l^(5)s( 3 - p i 11 )| 2 =- 327r3a3 -= [256m 8 (124s + 21*) - 64m 6 

' q(q)9K 1 ' 405m c 5 (4m2 - s) (4m2 - *) 5 L cV 7 

(288s 2 + 341s* + 63* 2 ) + 16m 4 (164s 3 + 456s 2 * + 321s* 2 + 63* 3 ) - 4m 2 *(106s 3 
+190s 2 * + 115s* 2 + 21* 3 ) + 11s* 2 (2s 2 + 2s* + * 2 ) 


q(q) + g ->• cc( 3 Pj 1] ) + q(q): 
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I ^ W 3 ^ 1 ')! 2 = 


327ra 


3^,3 


614400m( 2 - 34816m 40 (15s 


2025m®t (4 m 2 c — s) (4 rn? c — t) 5 L 
+13t) + 256m® (840s 2 + 1386st + 443t 2 ) - 64m® (420s 3 + 1708s 2 t + 1315st 2 + 177t 3 ) 

+16m 4 t (464s 3 + 1060s 2 t + 519st 2 + 43t 3 ) - 4m 2 t 2 (138s 3 + 206s 2 t + 87st 2 + 17t 3 ) 

+7st 3 (2s 2 + 2 st + t 2 ) 1 (41) 

q(<!) +9 cc( 3 Sf ] ) + q(q): 


|A r s (»)»( 3 Sf l )| 2 = 


7r 3 a 3 


. 128m® (20s 3 + 69s 2 t - 39st 2 + 20t 3 ) 

81 mist (4 m 2 c - s) (s + t) 3 L c v ’ 

-32m 4 (40s 4 + 113s 3 t + 27s 2 t 2 + 10st 3 + 20t 4 ) + 4m 2 (108s 5 + 193s 4 t + 41s 3 t 2 
+225sV + st 4 + 20t 5 ) - 11s (4s 5 + 3s 4 t + 7sV + 7s 2 t 3 + 3st 4 + 4t 5 ) 1 (42) 

q(q) + g -> cc( 4 4 81 ) +q(q)- 


|A', ( ,- )i ,( 1 4 8| )l 2 


107r 3 a 3 [m 2 (44s 3 + 92s 2 t — 4st 2 + 44t 3 ) — 5s (s 3 + s 2 t + st 2 + t 3 )| 
27 m 3 (s — 4mj?) (s + t ) 3 (—4m;? + s + t) 


(43) 


<?(e) + 9 -> cc( 3 pf 1 ) + <j(e): 




27r 3 a 3 


512m® (5s 2 + 26st + 25t 2 ) 


9m® (s — 4m;?) (s + t ) 4 (—4m); + s + t) 

+64m 4 (s 3 - 23s 2 t - 111st 2 - 19t 3 ) + 4m 2 (s + t) (57s 3 + 169s 2 t - 3st 2 + 61t 3 ) 

—31s(s + t ) 2 (s 2 + t 2 ) 

?(?) + 9 ->• cc( 3 SJ 81 , 3 d[ 8] ) + g(g): 


(44) 


|jv,(,-) 9 ( s s! 81 , 3 4 81 )l 2 = 


2Ti d a 


3^,3 


27-\/l5m®st(s +1 ) 4 L 
—32m 2 (5s 5 + 16s 4 t + 14s 3 t 2 + 14s 2 t 3 + 16st 4 + 5t 5 ) 

+5(s + t ) 2 (4s 4 - s 3 t + 8 s 2 t 2 - st 3 + 4t 4 ) 


128m 4 (5s 4 + 2s 3 t + 21s 2 t 2 + 2st 3 + 5t 4 ) 


(45) 


g + q ->• cc^P^) + g: 
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lA^f ? 1 )! 2 = 


1287 r 3 a 3 (12 m 2 c - s) (880m 4 - 112m 2 s + 13s 2 ) 
3645m 3 s (4m 2 — s ) 5 


x [2 1 (s — 4m 2 ) + (s — 4m 2 ) 2 + 2t 2 ] 
g + q —)■ cc( 3 P 1 [1] ) + g: 


( 46 ) 


| A '» S ( S P |' 1 )| 2 = 


2567T a 


3„,3 


64m®(31s + 84 1) - 16 m 4 (73s 2 + 150st 


1215m 3 (s — 4 ml) 5 
+841 2 ) + 4m 2 s (53s 2 + 88 st + 66 t 2 ) - 11s 2 (s 2 + 2 st + 2t 2 ) 

g + q —>■ cc( 3 P^) + g: 


(47) 


|JV »-»(+ 1| )| 2 = 


256vr 3 a 3 


92160m 40 - 512m 8 (71s + 90 1) + 64m® 


6075m 3 s (s — 4m 2 )" 

(33s 2 + 404st + 180t 2 ) - 112m 4 s (s 2 + 46st + 32t 2 ) + 4m 2 s 2 (33s 2 + 112st + 98t 2 ) 
-7s 3 (s 2 + 2 st + 2 1 2 ) 1 (48) 

g + g -+ cc^Sf 1 ) + g: 


|jv«( 3 s ! 81 )| 2 = 


87t 3 q; 


3^,3 


45056mi° 


243m 3 t (4m^ — s) 3 (—4m^ + s + t) L 
—256m 3 (84s + 205 1) + 64m® (128s 2 + 496st + 475t 2 ) - 16m 4 (224s 3 + 780s 2 t 

+1029st 2 + 540t 3 ) + 4m 2 (180s 4 + 676s 3 t + 1155s 2 t 2 + 936st 3 + 270t 4 ) 

-11s (4s 4 + 17s 3 t + 35s 2 t 2 + 36st 3 + 18t 4 ) 

q + q -+ cc^sf 1 ) + g: 

1 [s]v | 2 _ 4007t 3 q; 3 [16m 4 - 8m 2 (s + t) + s 2 + 2 st + 2t 2 ] 


(49) 


|tf«W J )r = 


81m 3 s (s — 4m 2 )^ 


(50) 


g + g —> cc( 3 Pj 8 ') + g: 


l JV 5a( 3p j l )| 2 = 


— 167T 3 CK 3 


29440m 8 - 128m® ( 68 s + 115t) + 32m 4 (59s 2 


27m 3 s (s — 4m 2 ) L 
+205st + 115t 2 ) - 8 m 2 s (64s 2 + 121st + 90t 2 ) + 31s 2 (s 2 + 2 st + 21 2 ) 


(51) 
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q + q —> cc^Sj 8 ', 3 Df ') + g: 


|Wm( 3 MVM 81 )| 2 =- 7 =-- 81920m, 

81\/l5 m\t (s — 4 ml) (—4 m 2 c + s + t) L 

—1024m 4 °(80s + 157 1) + 256m 3 (140s 2 + 569st + 535t 2 ) - 128m®(80s 3 + 425s 2 t 
+710st 2 + 378 1 3 ) + 32m 4 (70s 4 + 353s s t + 705sV + 630st 3 + 1897 4 ) - 4m 2 s(80s 4 
+353s s t + 700sV + 684st 3 + 252t 4 ) + 5s 2 (4s 4 + 17s 3 t + 35 s 2 t 2 + 36st 3 + 18t 4 ) 1 (52) 


167r 3 « 3 
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